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Abstract 

We introduce a large version of the spin quantum Hall transition 
problem. It is formulated as a problem of Dirac fermions coupled to 
disorder, whose Hamiltonian belong to the symmetry class C. The 
fermions carry spin degrees of freedom valued in the algebra sp{2N), 
the spin quantum Hall effect corresponding to = 1. Arguments 
based on renormalization group transformations as well as on a sigma 
model formulation, valid in the large N limit, indicate the existence 
of a crossover as N varies. Contrary to the N = 1 case, the large 
N models are shown to lead to localized states at zero energy. We 
also present a sigma model analysis for the system of Dirac fermions 
coupled to only sp{2N) random gauge potentials, which reproduces 
known exact results. 
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1 Introduction. 



The landscape of delocalization transitions is wider in two dimensions than 
in higher dimensions, parallehng the classification of the new random ensem- 
bles of ref.[0. Many of these transitions may be modeled by Dirac fermions 
coupled to various disordered variables. When formulated as field theories 
these transitions are usually mapped into difficult strong coupled systems re- 
quiring the identification of non-trivial infrared fixed points. In this respect, 
the recently introduced su{2) spin quantum Hall transition 0, which corre- 
sponds to random Hamiltonians belonging to class C of ref . , appears as an 
exception. Indeed, using a network formulation, it has been argued |0 that 
the critical properties of the latter model are potentially described by perco- 
lation in two dimensions. This opened the possibility of exact computations 
of characteristic critical exponents 0] or of the mean conductance 0. See 
also ref. for a spin chain formulation of this model. However, these results 
still resist to a field theory derivation, see eg. . Such alternative approach 
would be useful for the related, but yet unsolved, problem of the quantum 
Hall transition. 

The aim of this paper is to study a large version of the spin quantum 
Hall transition. We shall formulate it as Dirac fermions with sp{2N) spin 
degrees of freedom whose Hamiltonians belong to the class C. This require- 
ment, which forces us to choose the algebra sp{2N), is compatible with the 
introduction of four types of disordered potentials. These potentials are all 
generated under renormalization group transformations. Some of the results 
valid in the = 1 spin quantum Hall effect extend to the large N gen- 
eralization. This holds true for instance for the spin-charge separation for 
particular fine tuned versions of the models. However the following analysis 
indicates that there exists a crossover as N is increased. Namely, we will find 
that, contrary to the = 1 case, zero energy states are localized at large N. 

We analyse these models in two steps. First, renormalization group trans- 
formations allow us to identify the field theory describing the universality 
class of such models. It corresponds to disordered variables isotropically dis- 
tributed among the four possible types of disordered potentials. We then 
use supersymmetric sigma model techniques to analyse the universal model, 
which for class C Dirac fermions is a sigma model on the Riemannian sym- 
metric superspace OSp(2\2) / GL{1\1) , of type i5III|CI according to ref.[0. 
This sigma model, characteristic for the symmetry class C, also appeared 
in the context of disordered superconductors with spin rotation symmetry 
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and no time reversal invariance |^. The resulting sigma model turns out 
to be a massive theory. As a consequence the low energy limit is non criti- 
cal. This is in contrast with class D Dirac fermions which were described in 
refs. [plot |ll[] by a massless sigma model of type CI | DHL In both cases, the 



derivation of the sigma model is justified only at large N . 

For completeness, we have also included a sigma model description of 
Dirac fermions coupled to sp{2N) random gauge potentials. This corresponds 
to a particular fine tuning of the previously discussed models enlarging the 
symmetry of the Hamiltonians which then belongs to class CI. Surprisingly, 
large sigma model analysis reproduces the exact result that the critical 
theory is an osp{2\2)k=_2N WZW model. The fixed point of class CI, includ- 
ing the WZW term, was also found in []T^ [T^ using replica and in by 
direct means. 

The paper is organized as follows. The models as well as their super- 
symmetric formulations are introduced in Section 2 and 3. There, we also 
introduce the appropriate orthosymplectic transformations and the effective 
actions. Symmetries of the pure systems are identified in Section 4. The beta 
functions computed in Section 5 show that the model with isotropic random- 
ness is universal and attractive at large N. The sigma model formulations of 
the latter, which is valid in the large limit, is analysed in Section 6. We 
show that, at large scale, this model is driven to a strongly coupled massive 
phase, due to the absence of non-trivial topological 6'-term. This leads to our 
main result concerning localization of zero energy eigenstates. Spin-charge 
separation in the pure system and sigma model formalism for perturbations 
by random gauge potentials, corresponding to random Hamiltonians in class 
CI, are presented in Section 7 and 8. 

Acknowledgments: It is a pleasure to thank M. Zirnbauer for exchanges 
on related matter, for pointing to us that spin-charge separation should ex- 



tend to any A^ and for furnishing us ref.[15|, and M. Bauer and A. LeClair 



for numerous discussions and related collaborations. 
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2 The models. 



We consider 2N species of Dirac fermions, coupled by disorder in class C |T|. 
The Dirac Hamiltonian for generic disorder in this class is 

fa + M + m -2id + A\ 
\-2id + A a-M-m) 

where d = {dx—idy)/2, d = {dx + idy)/2, A = A^ — iAy, A = A^ + iAy. Here, 
Af^ = J2a ^fii^y y)^"' random sp{2N) gauge potentials, a = ^^(x, ?/)r" 
is a random "spin" potential, and M = Y2-Mi{x,y)T^ and m = m{x,y)l 
are random mass terms. We consider generators of sp{2N) defined by the 
relation 

where the superscript T denotes usual transposition, and S is the symplectic 
unit 

■ I \") ■ 

The generators T*, with the property T* = S(T*)"^S~^, form the complement 
of sp{2N) with respect to sl{2N). The generators and T*, together with 
the identity 1, span the algebra gl{2N). Remark that for = 1 there is no 
generator T\ so that sp{2) ^ su{2). 

The Hamiltonian (|I]) enjoys the following particle-hole symmetry defining 
symmetry class C (spin rotation invariance and no time reversal symmetry) 
in the classification of [|l| 

H = -CH'^ , (3) 

with C = o"i ® S an antisymmetric matrix. This relation implies that the 
eigenvalues occur in pairs with opposite signs, and it relates the advanced 
and retarded Green functions 

Gr{E) = -C GAi-Ef . 

We consider centered gaussian distributions for the four types of disorder, 
with strengths qa, Qa, Qm and (positive real numbers for real disorder). 
The disorder variances are: 

{m{x)m{y)) = ^ 6^'\x - y) , {M{x)M{y)) = ^ 6^'\x - y) 
{a\x)a\y)) = 5'^' 5^'\x - y) , {A'^ix)A\y)) = ^ 5'^'5^'\x - y) 
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The reason to consider the four types of disorder is that they are generated 
by renormahzation of the effective action. As we will see in the following 
sections, there are particular cases where the symmetry of the Hamiltonian 
([l|) is greater than (|^) and the renormahzation group flow closes on a subset 
of the disorder coupling constants; this is the case for example for the random 
vector potential alone, when the symmetry class changes to CI.. 

The single particle Green functions are defined by the functional integral 
J D\E'*D\E' exp(— S") with Z the partition function and 

S= f— (H - S) ^>{x) (4) 

J 27r 

where £ = E + ie. For e = 0~^, this defines the retarded Green function 

Gr(x,x';E) = lim (x\ -\x') = lim i(^(x)^*(a;')) (5) 

£^0+ H — [E + is) e^0+ 



Letting 



^ = ( ), vI/^ = (^-,^-) (6) 



where ip^ is a 2N-component fermion ip^^^i = 1, . . . ,2N, and similarly for 
one finds 

/(Px 

+ im{il)-i)+ - iJ-ip+) + i{tp-Mij+ - ij_Mip+) - iS^E^ (7) 
where = ip-^+ + 



3 Supersymmetric effective action. 

Since the disorder is gaussian and the deterministic part of the hamilto- 
nian is free, we can use the supersymmetric method to compute disorder 
averages of Green functions. For each fermion field, we introduce bosonic 
partners, i = 1, . . . , 2N, so that the inverse of the partition function 

for fermions become a partition function for bosons 

Z{a,m,M,A)-^ = [ Dp e-^^^^f^^ (8) 
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In order to simplify the notation and to take advantage of the symmetries of 
the problem, we are going to introduce supermultiplets, each containing AN 
fermion fields and 4N boson fields 

















= 











with S defined in eq. (^. Here the index i = 1, . . . , 2N was omitted, and 
the superscript T represents the usual transposition, tp^ is a column vector 
and tp- a row vector. We define also orthosymplectic transposes 0* of and 
0* of by 



The inner product associated to this transposition is skew 



(9) 



Similarly, 



(10) 



One can define a transpose of 8A^ x 8A^ supermatrices A by (A0)* = 0*A*. 
Its explicit relation to the usual supertranspose (we use same symbol T to 
denote the transposition for the usual matrices and the supertransposition 
for the supermatrices), is given by 



A' 
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S 



E 



FF 





-S 



E 



BB 



EpF and Ebb being the projectors on the fermion-fermion space and boson- 
boson space respectively. The present definition of the supermultiplets 
and 0*, as well as of the orthosymplectic transpose is different from the one 
used in |TT]], and it is adapted to the present type of disorder. To be able to 
compare the symmetry properties of supervectors and supermatrices in the 
two cases it useful to retain the following relations 



^class C 



Lsusy V-'class D 



class C 



class D 





s-1 



■-susy ) 



(12) 
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and similarly for and 0*. 

The free theory of fermions plus the bosonic ghosts is conformal, with 
Virasoro central charge c = and with action 

5'cft = 2 y ^ {i^-d,i^+ + i^-d,i^+ + P--d,P+ + P-dzP+) (13) 

^ (0* 9,0 + 0* dj) . 

The short distance singularity of the holomorphic fields in the supermultiplet 
can be written as 

(p(zW(w) - (14) 
z — w 

and similarly for the supermultiplet 0. 

After performing the integrals over the disorder, we obtain the following 
effective action 

1 f cPx 

SeE = Scft + 7777 / TT- i9a Oa + gmOm + Qm Om + QA O a) , (15) 

/iV J /vr 

with the following operators perturbing away from the conformal field theory 

=(0*r»0)(0*r''0), (16) 
Om =^(0*0)(0*0), 

Oa =(0*r''0)(0V'^0). 

It will turn out to be useful to know alternative ways of writing these oper- 
ators. Using the antisymmetry properties (|T^) as well as the cyclicity of the 
sup ert race, one has: 

C„ = STr( 00* 00* ) = STr( 00* 00* ) (17) 

Om = ^STr( 00* 00* ) = -^STr( 00* 00* ) 

Om = STr( 00* 00* ) = -STr( 00* T' 00* ) 

Oa = STr(00*r"00*r") . 
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Finally, one can consider the energy term, which can be written in the 
present notations as 

$^=(0*^3 0), (18) 

where S3 = I2N ® 0"3 ® Isusy In order for the integrals over the bosonic field 
to converge, one has to relate the bosonic fields to each other by complex 
conjugation 

pi=p+, pi=p-. (19) 

Then, a positive imaginary part of the energy insures convergence of the 
bosonic functional integrals. 



4 Symmetries. 

Before perturbation by the disorder operators (|T6|), the action ( [T3| ) possesses 
conformal invariance with current algebra symmetry [^. Bilinears in the 
chiral supermultiplets (00*) generate an osp{4N\4N) algebra. 

Two natural subalgebras of osp{4:N\4:N) are sp{2N) and osp{2\2). The 
sp{2N) algebra, which we refer to as the spin algebra, is generated by the 
currents, 

.r{z) = {(p'r''(P){z) (20) 

They satisfy 

fab 

r{z) J^(0) ~ ^ J"(0) + reg. 

z 

with /"^ the sp{2N) structure constants, showing that they form a represen- 
tation of the affine algebra sp{2N) at level zero. 

The osp{2\2) algebra, which we refer to as the charge algebra, is generated 
by the sp{2N) scalar 

K{z) = Tigi^2N) (0 0*) {.) = Trgi(2iv) {Ej (0 0*) (,) E^) (21) 

where the trace Trgi(2N) is over the spin indices and Ei is any orthonormalized 
basis of gl{2N). A convenient basis is be made of r", T* and the identity 
1. The currents K may be viewed as a 4 x 4 supermatrix, Ka/3 = J2i 'Pia^'lp-i 
with components: 

J = i^W. , J± = S,,V'iV'i (22) 
H = (3X131 , 5± = 7/>i/3^ , 5± = S,,/5i/3i 
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The even part of K is made of two blocks KpF and Kbb with: 



The currents (^) generate an osp(2|2)fc current algebra at level k = —2N, cf 
ref.0. In particular, the Kbb block generates a so(2) algebra and the Kpp 
block generates a sj9(2) algebra. 

The osp{2\2) generators commute with the sp{2N) currents: 

[r,K] = o. 

The supermultiplets (p and 0* transform as affine primary fields with value 
in the tensor product of the defining vector representation of sp{2N) by the 
4-dimensional representation of osp{2\2). One may alternatively |TB| view 
the field as a rectangular supermatrix with components (j)ia, i = I, - ■ ■ , 2N, 
a = 1, ■ ■ ■ , 4, on which the algebra sp{2N) acts by multiplication on the left, 
while osj9(2|2) acts by multiplication on the right. These actions of course 
commute. 



5 One loop beta functions. 

The first step in studying the effect of the disorder perturbation in (|l^) is to 
analyse the renormalisation group flow for the coupling constants. Even at 
one loop, beta functions code for important properties of the system. They 
can be deduced from the operator product expansions (OPE) of marginal 
perturbing operators Oi, where i stands for one of the indices of the disorder 
operators ([161) . Using the following notation for OPE 

zz 



the one loop beta functions are given by [T^ 

1 



(3k = ldigk = -^l^C.jgigj 



After some careful calculations, we obtain the following results 



{2N + 1)/ . . (N-l) \ 1 



= — — y^gAigm + go) - gmga - 



2 
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1/ ^2 1/2 2 N 2(A^+ 1) 2 1 , 

/3a = + +^(^a-£/A/) + 9a + j;p9a{9m- 9m) 

Pa = '29A{9a + 9m) + J^9a{9M - 9a) + ^(^9a + 2gA^i9m- 9m) 

2(A^ + 1) , , 1 , , 1 , , , , 

Pm = 9A{9a + 9m) + j;^9m{9m - 9a) + j;^9m{9m - 9rn) i^S) 

Readers who are interested in details can find some hints in the Appendices. 

From now on, we shall be interested in the large N limit. The beta 
functions ( p3D then simplify dramatically 

Pm = -'^9m{9a + 9m) + ^9A{9a + 9m) , pA = ^{9a + 9m)^ + '^9^ 

Pa = '^9A{9a + 9m) , pM = '^9A{9a + 9m) (24) 

Of course the large N limit only applies for g <^ N. 

After a large number of RG iterations, RG trajectories are asymptotic 
to directions which are preserved by the RG flow and which pass through 
the origin. One can show that there exist six stable directions. But only 
one is attractive in the region where coupling constants are positive. This 
direction, which we call the "isotropic direction", corresponds to the case 
where all coupling constants are equal: Qm = 9a = 9m = 9a = 9- The 
associated beta function is 

Pg = ^9' (25) 

To show the isotropic direction is attractive, we use the same method as 
in p. We project RG flow on the sphere and parameterize coupling constants 
with a radial coordinate p and three angle coordinates 6i. We rewrite the 
RG equations as 

e, = pp,{eue2,e,) and p = p'Ppie,, 62,63) 

We develop the Pi around the isotropic direction, and compute the corre- 
sponding eigenvalues. These eigenvalues are all negative, proving that the 
direction is attractive. Furthermore, it can be shown that p > 0. These 
results are corroborated by the all order computations of Appendix C. 

Thus, the isotropic direction describes the universality class in the region 
where coupling constants are all positive. It corresponds to a strong coupled 
system that we will study in details in the following section. 
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6 The sigma model approach. 



When the disorder couphngs are all positive (and is large), the renormal- 
ization group flow is attracted towards the line Qa = Qm = Qm = Qa = 9 
("isotropic disorder"). On this line, preserved by the flow, the low energy 
physics can be described by a non-linear sigma model with a topological term. 
In this section, we derive the effective action for this sigma model. Its target 
space is the Riemannian symmetric superspace 05'p(2|2)/GL(l|l), denoted 
Z}III|CI in 1^ (note that there exists another 05'p(2|2)/G'L(l|l) symmetric 
superspace, denoted CI | Dili, appearing in the context of the symmetry class 



Dm)- 

When all the couplings are equal, the perturbation term in the La- 
grangian, quartic in the Dirac flelds can be decoupled by a Hubbard-Stratonovich 
transformation involving an unique supermatrix fleld Q. Using the fact that 
(0* 0) = (0* T* 0) = and denoting 

5 = 00* + 00* 

the perturbation term becomes 

^pert = ^ ((0V«0)(0*r'^0) + ^(0*0)(0*0) (26) 

+ (0*T*0)(0*r^0) + (0*r'^0)(0*r'^0)) 

= — STr 5 f T^B + —B + rB T 
AN \ 2N 

The supermatrix B obeys the relation B = — -B*, relation which deflnes an 
element of the algebra osp{4N\4N). Conjugation by the generators of the 
subalgebra gl{2N) projects B on the identity on this subalgebra 

r'^B + ^B + r'B r = (Tr,,(2iv)5) ® I27V , 

so that the perturbation term reads 

Vrt = ^ STr (Tr,,(2^) (00* + 00*))' . 

This interaction can be decoupled using a supermatrix belonging to the 
osp{2\2) algebra, Q ~ Tigi{2N) (00* + 00*). The resulting effective lagrangian 
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IS 

Leff = 0*90 + 0*90 + STr {QTrgi^2N) (00* + 00*)) -jSTiQ 
it^t.iQ d\f<p\ N^^^. 



2 



a QJUJ (27) 

The next step is to perform the gaussian integrals over the Dirac fields, 
resulting in the following effective action 



S[Q] = -j j ^STrg2 + ArSTrln(^ 



where STr combines the operations of taking the supertrace STr over the 
matrix indices and integrating over position space. The factors S3 under 
the logarithm appear after the transformation 1— ^30 and 0* ^ 0*^3, 
correcting for the fact that the complex conjugate of is not 0* but 0*^3, 
see eg. ([T9|) . The number appears now as a factor in the action, suggesting 
to treat the integral in the saddle point approximation. The saddle point 
equation for the action ([28| ) is given by 



'^^"^-Tr,,,)(x|p-^|x> , with 2)=(^ ^) (29) 
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We look for a spatially homogeneous solution of the form Q{x) = nT,^. The 
saddle-point equation then reduces to 

Cutting off the integral in the ultraviolet by \k\ < I/Eq yields the equation 

l/2g = In (1 + (2/i£o)"') ~ -2 ln(2/i£o) , (30) 
and by inversion, 

/i = (2£o)-VVeV2. - 1 ^ (24)-^ e-'/'^ , (31) 

As the dynamically generated mass yU is a renormalization group invariant, 
eq.(pID corresponds to a beta function Pg = Ag^ in agreement with eq. (pSf) . 
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Symmetry of Q and the saddle point manifold. The bilinear in the Dirac 
fields B = (j)(f)^ + 00* belongs to the complex algebra osp(4A^|4A^) defined by 
the relation B = —5*. The even part of this algebra is sp{AN) ^ so{4N) , with 
so{4N) in the fermion-fermion {FF) block and sp{AN) in the boson-boson 
(BB) block (this is due to the fact that the matrix 7 in ([TlD is symmetric 
in the FF sector and antisymmetric in the BB sector). Taking the trace 
over gl{2N) in B leaves us with an object belonging to the osp{2\2) algebra, 
but now with the so(2) part in the BB sector and the sp{2) part in the FF 
sector. To see the way it happens, let us look at the FF part in a osp(4A^|4iV) 
matrix, defined by M* = —M. It is easily verified that this block has the 
structure 

where Mn, M12 and M21 are 2N x 2N ordinary matrices obeying M^^ = 
-M12 and = -M21. This implies that MpF belongs to so{AN). Let us 
now take the trace over the colour indices 

Mff = Tr,,(2;v)M^p = ^ . (33) 

Using the cyclicity of the trace and the antisymmetry of the matrix E, we can 
show that = M12 and = M21, which means that Mpp belongs to the 
sp{2) algebra. Similarly, Mbb can be shown to belong to so{2). The matrix 
Q inherits this symmetry from the object to which it couples, B = TTgi(^2N)B, 
therefore it belongs to osp{2\2). 

When decoupling the interaction part with the help of the supermatrix 
Q, one of the question that has to be addressed is the choice of the contour 
of integration. In particular, solving this question allows to choose the ac- 
ceptable solutions for the saddle point equation, that is the ones which lie 
on the contour of integration or which can be attained from it by analytical 
continuation. These questions have been addressed in detail in 0] for the 
class C. There, it was shown that the dominant diagonal saddle point is of 
the form 

Qo = fi Tj^ , 

with II3 = cr3(S)lsusy being an element of osp{2\2). Due to the global 05^(212) 
symmetry of the effective action, this saddle point extends to a saddle point 
manifold 
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where T is a constant element of OSp{2\2). Since the stabihzer of S3 is 
GL(1|1), the saddle point manifold is the coset space OSp{2\2)/GL{l\l). 
This coset space can be parameterized by T = expX, with {X, E3} = 0. 

The convergence conditions for the integrals over Q restrict the saddle 
point manifold to a real submanifold of the complex space OSp{2\2)/GL{l\l). 
In the FF sector, convergence of the integrals over Q can be insured by 
choosing Q^pp = Qff- At the level of the saddle point, this translates into 
Xpp = —XpF- Therefore, the fermion-fermion sector of the saddle point 
manifold is isomorphic to the compact symmetric space Sp{2)/U{l). The 
convergence conditions on Q in the boson-boson sector are more involved ; 
on the saddle point manifold they can be reduced to X^^^ = Xbb, showing 
that the bosonic part of the saddle point manifold is non-compact. When 
averages of n Green functions are considered, it is isomorphic to the non- 
compact symmetric space S0*{2n)/U{n), with S0*{2n) some real form of 
S0{2n). For n = 1, the boson-boson sector is empty. The two symmetric 
spaces form the base manifold of a Riemannian symmetric superspace of type 
Dm\Cl §. 

Gradient expansion. The next step in the derivation of an effective ac- 
tion is to perform a gradient expansion of the action p^). The low energy 
configurations are given by the slowly varying field 

q{x) = Q{x)/fi = T{x)J:,T{x)-' , (34) 

where T{x) is a (slowly varying) element of 05^(212). Note that q{x) satis- 
fies the nonlinear constraint g(a;)^ = 1. The degrees of freedom q{x) corre- 
spond to the Goldstone modes of the broken symmetry OSp{2\2) G'L(lll). 
Fluctuations transverse to the saddle point manifold are massive and can be 
neglected at this stage. 

The effective action for the Goldstone modes is a non-linear sigma model 
on the symmetric superspace OSp{2\2)/GL{l\l) described previously. This 
sigma model may support a topological term, since Il2{Sp{2)/U{l)) = Z. 
The easiest way to extract the coupling constants of the kinetic and topolog- 
ical term is by using the non-abelian bosonisation . In the supersymmetric 



setting, this method was used and explained in detail for class D in [|rT| and 



it can be applied with minimal changes to the present case. We want to eval- 



uate the action (|2q) on configurations of the type (|3J). Due to the nonlinear 



constraint g(x)^ = 1, the first term in ( pH]) vanishes. The second term can 



14 



be written, by undoing the integral over the Dirac fields, 



e-^t"?] = jv(l)V(p exp - J dPx {(p^dcp + 0*90 + /i0*g0 + ficp^qcp) . (35) 

Here q has to be understood as acting like the identity on the spin indices 
z = 1, ■ ■ ■ , 2N. The free Dirac theory plus the bosonic ghosts is equivalent to 
a WZW model with action 

Wospimm) [M] = ^ [ d^x STr (M-'d^My + , (36) 

where the matrix M takes values in a subspace of the complex supergroup 
0Sp{AN\4N), and the topological term is expressed by assuming some ex- 
tension M of M to a 3-ball B that has position space S for its boundary 
{dB = S) 

T[M]= [ d^xe^^xSTiM-^d^MM-^d^MM-^dxM . (37) 
Jb 

The rules of bosonisation for the last two terms in the exponent in (RH) 
are to replace the bilinears 00*^3 and ^300* by i~^M resp. i~^M~^, where 
the factor is a large mass scale, of the order of £q^, which enters for 
dimensional reasons. Up to a conjugation with the matrix diag(l, S), these 



are the same bosonisation rules as in |11|. The term STr [MTi^q + qH^M ^ 



can be viewed as a kind of mass term. At large it forces the field M to 
follow qT.3. This approximation is valid at momentum scale k -C (fx/iY^'^. 
Neglecting the fiuctuations we can set MTi^q = 1, which yields 

S[q] =2iViy[gS3]. 

where the factor 2N appears from taking the trace over the spin indices. 
Here WlqJ:^] is the WZW action on OSp{2\2). Recall that S3 G osp{2\2) 
so that ±zS3 = exp(±z7rS3/2) and gSs belong to OSp{2\2). However, gEa 
does not explore all this group as q is only a function on the coset space 
OSp{2\2)/GL{l\l). Evaluating the topological term for this configuration 
can be done by making a smooth extension of M = gEa to the ball B with 
radial coordinate < s < 1, for example 

M(x, s) = T{x) exp (±zs7rE3/2) T{x)-\Ti^3) . 
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At s = 1 we have M{x, 1) = g(x)S3, while for s = we get M{x, 0) = =F^S3, 
independent of x. Inserting this extension into the expression (|37|) for r[M], 
and converting the integral over B into an integral over S = dB, we find a 
theta term 

i If 

— r[gS3] = ±— J d?xe^^ STrqd^qd^q = ±Stop(g) • 

Since the value of the WZW topological term does not depend on the ex- 
tension, the two opposite expressions for ir[gS3]/247r are equivalent and 
Stop{(l) G ivrZ. Gathering the kinetic and topological term, we obtain the 
following effective action 

2N f 2N f 

S[q] = -TF~ / d^x'^Tid^qd^q±— / (fxe^u^T^qd^qduq. 

g=rs3T-i 167r J 32 J 

(38) 

The angle of the theta term is 6' = ±2A^7r. It contributes trivially to the path 
integral as the topological action is multiplied by 2N so that 2A^5'top G 2i7rZ. 
The effective action is thus: 

2N f 

Ses[q] = -— I d\STrd,qd,q . (39) 

The natural ultraviolet cut-off for this effective action is ~ 2£oe^/^^ since 
in deriving it we neglected transverse modes of effective mass /x. 

Action (|39D is conjectured to be a massive theory as it is a sigma model 



on a symmetric space with positive curvature. Recall [T^ the one loop renor- 
malization group equations for sigma model metrics Gab'- 

Idl Gab = ~Rab 

with Rab the Ricci curvature. For symmetric spaces the Ricci tensor is pro- 
portional to the metric. We need to compute this proportionality coefficient 
in our case. Since q = TT^^T'^, the tangent space at the point g = E3 
is spanned by elements of the form [X, S3] with X G osp{2\2) / gl{l\l) . By 
construction we can choose X such that {X, S3} = 0. The metric is then: 

GiX,X) = -^STr([X,S3])2 = ^STr(X2) 

where the supertrace is understood in the defining 4 dimensional represen- 
tation, and A = 1/2N. Similarly, the Ricci tensor at g = S3 is defined by 



R{X,X) = -STr(adX)^ 
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Here the supertrace is in the adjoint representation. Recall that Ricci ten- 
sor are invariant under metric dilatations. To compute the proportionality 
coefficient we pick a particular element of osp{2\2) anticommuting with S3, 
e.g. X = EpF ® CTi- We have STr(X^) = —2. Diagonalizing the adjoint 
action, we get a set of eigenvalues zero, two bosonic non-degenerate eigen- 
values ±2 and two fermionic eigenvalues ±1 with multiplicity two. Hence, 
STr(adX)2 = 2(2^ - 2) = 4. Thus G = R/AnX. The RG equation then 
becomes: 

Wi A = +47rA^ 

At large distance, the model is driven to strong coupling, and it presumably 
becomes massive because there is no contribution from the topological term. 

The generated mass scale is of order ttitv — fie"^^^'^ <^ /i as the coupling 
constant is equal to at the ultraviolet cut-off fi. This is the energy scale 
at which the coupling constant A becomes of order one. 

As a consequence, the infrared fixed point is trivial and the zero energy 
states are localized with localization length of order l/m^r- In the localized 
regime, the behavior of the density of states is expected to be governed by 
the class C matrix ensemble 



7 Spin-charge separation. 

The conformal field theory with action ([131) admits a spin-charge separation. 
Its stress tensor can be decomposed into the sum of the Sugawara stress 
tensors associated to the spin and charge current algebras: 

Tcft = Tsp(2N)o + Tosp{2\2)_2N (40) 

Both Sugawara stress tensors have Virasoro central charge zero and are bi- 
linear in the corresponding currents: 

Tospi2\2)^2N = ~ 8(iV+ 1) 

The normalization of Tosp(2\2)^2N found in ref. [^. Eq. (^0|) is proved 

in Appendix B. 
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This may be checked by computing the dimension of the supermultiplet 
0. In the free theory, its conformal dimension is 1/2. The dimensions in the 
spin sector are Asp(2N)o = 2IN+1) '^i^h Cas. the casimir of the corresponding 
representation of sp{2N), cf. ref. I^Sj. For the vector representation this 
gives Asp(^2N)o = ^F+iy- charge sector, regular representations of 

osp{2\2) are labeled by two integers j, b and their conformal dimensions are 
^osp(2|2)_2Jv ~ '^2{N+i) ' [0- For the 4-dimensional representation with 

j = 1/2, b = this gives Aosp(2\2)_2N — I(W+T)' should, the spin and 

charge conformal dimensions add up to 1/2. 

For = 1, it was shown in ref. ||2^ that the four point correlation func- 
tion of the supermultiplet (j) may be factorized as the product of correla- 
tion functions in the sp{2N)q and osp{2\2)_2N conformal theories. However, 
this spin-charge factorization possesses peculiar properties inherited from 
indecomposability properties of representations of osp{2\2). In particular, 
osj9(4A^|4A^) decomposes as: 

osp{4:N\4:N) = osp{2\2) ® [1] + [8] ® sp{2N) + [8] ® [R] 

with [R] a {2N + 1)(A^ — 1) dimensional representation of sp{2N) and [8] 
isomorphic to the adjoint representation of osp{2\2). The spin currents J° = 
(j)^T°'(f) belong to [8] ® sp{2N) with [8] an eight dimensional indecomposable 
representation of osp{2\2). Thus, although these currents commute with the 
osp{2\2) charge generators they do not belong to a trivial representation of 
the charge algebra. 

This separation between spin and charge degree of freedoms still holds 
in the perturbed theory provided one fine tunes the coupling constants such 
that ga + dm = da + Qm = 0. In this case: 



-^charge = 9a {C^a " C^m ~ C^Al) 

= STr [ ( Ty<j>'T^ + ^00* + r>0*T^ ) 00* ] 

where we used again the antisymmetry property ([10|) . The sum in the above 
r.h.s. projects out the sp{2N) colour indices leaving only the osp(2|2)_2Ar 
currents (0). Thus 

Wgc = ^7.STr(iri?) (41) 
The remaining perturbing operator, 

L,pin=9AOA=gArr (42) 
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describes a current interaction involving only the sp{2N)q generators. Hence 
the total perturbation, 

Lpevt -^spin ~l~ -^charge 

is then the sum of the two commuting current-current interactions. 

The spin-charge separation can also be seen on the beta functions. When 
fine tuning coupling constants such that Qa + Qm = Qa + Qm = 0, the beta 
functions ( p3D decouple: (3a = '^^^n^^ 9a Pa = —(3m = —(3m = —jf9a- 
the fine tuned regime, disorder in the gauge potential A is marginaly relevant 
while disorder in the spin potential a is marginaly irrelevant. 



8 Sigma model approach for the spin-charge 
separated system. 

The line Qa = dm = Qm = 0, Qa = 9 > is stable for the renormalization 
group flow. It is attractive in the fine tuned regime ga + Qm = Qa + Qm = 0, 
ga > 0, and along it the flow is towards strong coupling, qa ^ oo. This 
model was formulated and analysed in [|12| using replica and in [jl4| by direct 



means or with supersymmetry. Comparison of the various methods was done 
in 



25|. Based on the spin-charge separation, it can be deduced that the 
low energy physics on this line is given by an osp{2\2)k=_2N theory. Let us 
derive the same result by using the sigma model approach. We shall obtain 
that the effective action is a sigma model on the supergroup OSp{2\2), en- 
dowed with a WZW term, the coupling constants being such that k = —2N. 



Strictly speaking |ri| , the resulting WZW model is defined on a submanifold 
of the complex supergroup which is a Riemannian symmetric superspace of 
type D\C (meaning that the bosons have an orthogonal structure and the 
fermions are symplectic). Let us remind that by bosonising the free Dirac 
fermions/bosons, one obtains a WZW model on a Riemannian symmetric 
superspace of type C\D, at level k = 1. Since the metric changes sign when 
passing from a space of type C\D to one of type D\C, the WZW model on 
the space of type D|C is well defined for negative values of the level k. 

On the fixed line ga = gm = gu = 0, (^^ 7^ the symmetry is not that of 
class C any more but that of class CI. The reason is that the Hamiltonian 
has now an extra symmetry 

H = TH^r-^ with r = i(T2 ® S , 
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which can be interpreted as a time reversal symmetry. 

The disorder perturbation is now simply the sp(2iV)o current-current per- 
turbation. Adding formally terms which are zero, we obtain 

Lspin = |^(0V'^0)(0V'^0) = |^STr#*E,(00*)E, (43) 

= 1^ STr (Trgi(2iv)00*)(Trgi(2iv)00*) , 

with Ej the (orthonormal) generators oigl{2N). Since {(fxp^Y = —(jxf)^, we can 
decouple this interaction by introducing a supermatrix Q ~ 1i gi(2N)<t'4'* ■ We 
can define an orthosymplectic transposition for Q hj = Trg;(2Ar)(l ® Q f. 
Remark that in contrast to the preceding section Q has no specific symmetry 
properties. In particular Q ^ — Q*, so the supermatrix Q belongs to a space 
larger than osp(2|2). After decoupling of the interaction term, the effective 
lagrangian becomes 

2N 

Leff = 0*90 + 0*90 + 2STr (gTrgi(2iv) (00*)) + STr QQ* 

9a 



Here, we have embedded Q in the osp(4|4) algebra represented by matrices 
of the form 

i = ^) ' ^^^^ A' = -A, D' = -D . 

It is easy to see that the diagonal blocks of A span two commuting osp(2\2) 
algebras, so we conclude that the space to which Q belongs is isomorphic to 
the complement of osp{2\2) © osp{2\2) in osp(4|4). 

In the absence of the energy term, the lagrangian (|4^) is invariant under 
the holomorphic/antiholomorphic transformations 

<P-^gLiz)<P , 0*^0*^^i(^), (45) 
4>^gR{z)(f) , 0*^0*^^^(^), 

Q ^ 9l{z) Q g^\z) , ^ g^(z) g-\z) , 

with giiz) and gR{z) elements of the 05*^(212) group: g\^ = gj^^ and (7^ = 
g^^. It defines an action of the group G = OSp{2\2) OSp{2\2) on Q. We 
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therefore expect conformal invariance of the osp{2\2) (or charge) part of the 
theory. This is not surprising, in view of the spin-charge separation, since 
the charge sector is unperturbed by the disorder. 

Integrating out the Dirac fields, we obtain the effective action 



Again, for large, the integral over the Q field can be treated in the saddle 
point approximation. The diagonal saddle point solution is also given by a 
constant matrix of the form Q = fiA^3- The saddle point equation fixes the 
value of fiA, 



to the fact that, at large A^, (3a = '^q'ai /^a is also invariant under the RG 
flow and it has the properties of a dynamically generated mass. 

As usually, the saddle point solution extends to a saddle point manifold, 
due to the invariance of the effective action. This time, the saddle point 
manifold is generated by 



which is a much larger manifold than that of constant matrices. Q satisfies 
the non linear constraint QQ* = —1- Remember that S3 = —S3. The 
stabihzer H of S3 is constituted of elements dia.g{gL, gn) obeying 



This is compatible with the multiplication law as ghhh = S3 gnhji S3 
for any pairs {gL,gR) and {hi^hR). This means that H ~ 05*^(212). The 
fluctuations around the saddle point are described by a sigma model on 
the coset space G/H = OSp{2\2) ® OSp{2\2) /OSp{2\2). Elements of G/H 
are pairs {gLi9R) with the identification {gL^gn) ~ {Qi^zh^z ,9Rh) with 
h G OSp{2\2). Convergence of the integrals on the saddle point manifold is 
insured, in the FF sector, by Q^p = —Qpp- In the BB sector, convergence 
requires g[ j^ = T^^g^^^^T.^, or Qbb = gL,Bg\,B^^- This is equivalent to 
demand that the hermiticity condition (|l^) on the bosonic component, 0^ = 
(^^S3, is preserved by the transformation (^Sf). The even part (or the base) of 
the saddle point manifold is then equivalent to the product of two symmetric 





Q = fJ'A 9l{.z) S3 g^i^z) , 



Ql S3 = S3 =^ gL = S3 gji S3 . 
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spaces, the fermionic sector being compact and the bosonic one being non- 
compact. The full saddle point manifold is a subspace of the complex group 
05'j9(2|2), being a Riemannian symmetric superspace of type D\C. 

In order to perform the gradient expansion, we are making use again of 
the non-abelian bosonisation. First, we redefine the supermatrix field Q{x) 

Q{x) = fiA gix)!:^ , with g = gJl^gRY.^ . 

The nonlinear constraint on the new field is g^{x) = g{x)^^ which is nothing 
else that the defining relation for an element of the (complex) supergroup 
OSp{2\2). This shows that the coset space G/H is diffeomorphic to OSp{2\2). 



When inserting this expression in the effective lagrangian (|4) the last term 
vanishes. The other terms can be bosonized using the same rules as in the 
preceding section. The free Dirac part gives the WZW action (p6D, while the 
"mass" term can be bosonised replacing again (fxp^T,^ and by i~^M 

resp. i~^M~^ = i~^M^. The effective action becomes 

^STr {g + M g-') . 

The mass term forces M to follow g. Neglecting the fluctuations of M around 
the minimum M = g, we obtain 

S[Q] = 2N Wosp(2\2)[g] , (47) 

where the factor 2N appears from taking the trace over the spin indices, on 
which g acts as the identity. The topological WZW term survives to this 
reduction from OSp{4:N\4:N) to OSp{2\2). The reality conditions discussed 
above mean that the restrictions gp and gB of g in the FF and BB sectors 
satisfy gp^ = gp gs = gs' respectively. This ensures the stability of the 
action (^). 

As already mentioned, Wosp{4:N\4:N)[M] is defined on a Riemannian sym- 
metric space of type C\D, while W^osp(2|2)[5'] is defined on a Riemannian sym- 
metric space of type D\C. The quadratic form defining the metric changes 
sign between the two types of spaces, since STr = Tybb — T^t^ff = Trgp— Tr^o 
in the first case and STr = Tr 50 ~ Tr5p in the second case. In order to have a 
well defined functional integral, the level k of the WZW action has to change 
sign between the two cases. We conclude that the action (|^) corresponds to 
an osp(2|2)fc=_2Ar theory. It is surprising that the saddle point approximation 
is able to reproduce this exact result. 
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9 Appendix A: Algebraic coefficients 



The aim of this appendix is to define the algebraic conventions used for 
the calculations. We will also gather various algebraic formulas helpful for 
reproducing our results for the beta functions. 

The ensemble of sl{2N) generators are denoted by T"^, and they can be 
divided in generators of the subalgebra sp{2N), r", and generators of the 
complement of sp{2N) in sl{2N), T\ The normalization we choose is 

The various structure constants are defined by 

where [., .] stands for commutators and {., .} anticommutators. 

When evaluating beta-functions, a certain number of algebraic identities 
are needed to simplify expressions. Here is a fist of algebraic coefficients that 
are needed : 



Coefficient 


Identity 


Value 


C 


rparpA ^ Q 


(4Ar2 - l)/2N 


C 


r^r'' = C 


(2iV + l)/2 


c 


T'T' = C 


(2Ar+ l)(iV- 1)/27V 


XtTT 




-l/2Ar 


Xttj- 




-1/2 


XTtT 


rvr* = XrrrT" 


{N -l)/2N 


XtT t 


t'^TW = Xrrr'^' 


1/2 


xrrr 


T'T^T' = xrrrT^ 


-{N + l)/2N 


Dj"]^ 


fCDAfUDB ^_jj^^^AB 






jcdajcdb ^ -Drr^"" 


2(iV+l) 


Drr 


jijajijb _ —£)^^§ab 


2(A^- 1) 


Drr 


jaki jakj _ —JJ^^^V 


2N 


Brr 


fjabifjabj ^ -BrrS'^ 


-2{N + l) 


Brr 


^cia^cib ^ -BrrS"" 


-2{N^ - 1)/N 


X 




{2N+1)/N 
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10 Appendix B: Stress tensor factorization. 

The decomposition may be verified directly by computing the operator 
products defining the stress tensors using Wick theorem with the normahza- 
tion ([T^). For the sp{2N) sector we get: 



4(A^ + 1)''^ 8(A^ + 1) 

where we used 0*0 = 0, by antisymmetry, and the value of the sp{2N) 
Casimir in the defining representation, r'^r" = (2A^ + l))/2. Similarly for 
the osp{2\2) sector we first introduce the appropriate basis of gl{2N) to 
decompose Tosp(2\2)^2N follows: 

Tosp{2\2),2N = - g^^"*"^ ^-^ lim STr ((00*)(^) Ei {(pcj)*) i^u,)E i) 

Because of the antisymmetry property, 0*0 = (j)^T^(f) = 0, eq. (p!0D, the two 
last terms are proportional to (0*90). The first term gives a contribution 
proportional to (0*90) but also to (0*r'^0)^. Gathering these contributions 
and using again the value of the sp{2N) Casimir as well as T*T* = {2N + 
1){N-1)/2N, we get: 

Tospi2i2,-.. = -4(]vTI) - ^j^)^^'-'^^' 

Hence: 

Tsp{2N)o + Tosp{2\2)_2N = -2{4>^9(j)) = Tcft 

which proves the claimed factorization. 
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11 Appendix C: All order beta functions. 



In this appendix we derive expressions for the beta functions using formula 



suggested in ref. p6|. These formula apply to current-current perturbations 



of WZW models of the form: 

where J" and J"^ are the left and right conserved currents of the WZW models 
generating an affine (super) algebra at some level k. In the case of a super- 
algebra the currents J" possess a bosonic or fermionic character depending 
whether their degree \a\ are zero or one. In our case, the underlying algebra 
is the affine osp{AN\4N) at level k = 1. The currents are bilinear in the 
supermultiplet J" = with X"" generators of osp{4N\4N). 

For the theory to be perturbatively renormalizable, one needs to choose 
the tensors d!f^ such that: 



■ \^K jL facfbd /^KLjI 

^ab^cdJi Jj — ^ij 

V (^ai^^bj — ^ba 



where rj"'^ is the Killing invariant bilinear form of the superalgebra and /"^ 
its structure constants. These conditions are satisfied by the four operators 
Oa, Oa, Ora and Om defined in eq. (|T^). 

With an appropriate renormalization prescription, the proposed beta 
functions [pHl read for = 1, 



= -C{h!, h!){\ + D{}if) + 2C{h'D{h), h'D{h)) - 2h'D{h)RD{h) (48) 

where D{h) is the matrix D{h)^ = D^^hj, C{x,y) is the row vector with 
components C{x,y)K = C^^xlVi, and finally 



h' = h{l - D{h) 



2\-l 



We shall assume that these formula are correct and capture perturbative 
contributions to all orders. 

With the normalization of eq. ([T3|) , = Qk/'^N . 
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One has to compute all the tensors Cf^, Df^ and i?f . The first one 
codes the operator product expansion of the perturbing operators: 

C^(z)C^(0) ^ 5Z Cf-O^^)) + reg 

c 



The coefficients Cf have been determined when computing the one- loop 
beta functions, see eq. fpSl) . 

The tensor Df^ may be computed using a variation on the free field 
representation of the currents. Namely introduce auxiliary copies 0q,, a = 
0, 1 or 2 of the supermultiplets with two point function {(j)a<Pa)a = 1- Here 
(■ ■ ■)a denotes the expectation value in the auxiliary Fock space associated 
the the auxiliary field (pa- Let 

Then one has: 

Hence, Df^ is computable simply using Wick theorem. 

The last tensor Rf^ is computable by looking at the following operator 
product expansion: 

T^(z)C^(0) ~ ^(2Df ^ + MDf^)C^(O) + reg 



with T^(z) = d^^JHz)J' 



z 



The result is summarized in the following tables: 

L)-™ = D^;^ = = D^"" = 1/N, D^"" = {2N + 1){N - l)/N 
D""/ = 1, = = N + 1, Df/^ = (A^ + 1)(A^ - 2)/N 

= 2N + 1, D^'^^ = -AT + 1, 
j^MA = (AT - 1) /AT, D*^" = -AT, " = {N^ - 1) /N 

The non-vanishing ^ = 2Df^ + R§Df^ are: 

F^^ = 2Ftr = 2F//^ = 2F^r = HN + 1) 
F^"' = = FT = -2Fr = 4(2iV + 1) 

jpAm jpAm in Tpam Tpam In rpmA jprnm jpma Tp2M n / AT 

^A = = ^A = 1^ = "^A = '^m = a = "^M = ^/iV 

pAa ^ paa ^ 2, F*^^ = F*^" = F^'" = F§^^ = -2{2N + 1){N - 1) /N 
F^'' = AN, F^^^ = Ff^ = 2{N - 1) /N, F^^^ = Ff^ = A{N^ - 1) /N 
= -2{2N - 1), F^ = -2{AN + 3), F"^^ = -2{AN + 1) 
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For arbitrary A^, the beta functoins are then derived from eg. . In the 
N = oo hmit, they reduce to: 



Pm = Pa = '2{ga + gM)gA (49) 

{.ga + gM){,ga- gu f 



Pm = 4:gAiga + gm)- '^gmi.ga + gu) 



{ga - gu - 2) 



These are easy to integrate. In particular {ga — gM) is a RG invariant in the 
large limit while the beta functions for 2gA ± {ga + gm) are separated and 
quadratic. 

It is then simple to verify that the isotropic line ga = gm = gu = gA = g 
stable and attractive to all orders and f3g = Ag'^. The fact that it is quadratic 
is in agreement with eq. (|30| ) and it provides a tiny check of the all order beta 
functions. The isotropic coupling grows indefinitely with the scale. Of course 
the large N approximation remains valid only for g <^ N. It is possible to 
verify that all RG trajectories starting in the domain of positive couplings 
sufficiently close to the origin are asymptotic to the isotropic line at large 
distances. This confirms the one loop computation. 
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